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Abstract 

We describe the real quasi-exactly solvable spectral locus of the 
PT-symmetric quartic using the Nevanlinna parametrization. 
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Following Bender and Boettcher [3J, we consider the eigenvalue problem 
in the complex plane 

w" + (C 4 + 2b( 2 + 2tJ( + X)w = 0, w{te- wi/2±wi/3 ) 0, t +oo, (1) 

where J is a positive integer. This problem is quasi-exactly solvable [3J: there 
exist J elementary eigenfunctions w = p n (C) exp(— i( 3 /3 — ibQ, where p n is 
a polynomial of degree n = J — 1. 

When b is real, the problem is PT-symmetric. By the change of the 
independent variable z = i(, ([I]) is equivalent to 

- y" + (z 4 - 2bz 2 + 2 Jz)y = Xy, y(te ±ni/3 ) 0, t -> +oo. (2) 

Polynomial h in the exponent of an elementary eigenfunction y(z) is h(z) = 
z 3 /3 — bz. The spectral locus Zj is defined as 

{(6, A) E C 2 : 3 y ^ satisfying ©}. 
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The real spectral locus Zj(R) is Zj D R 2 . The quasi-exactly solvable spectral 
locus is the set of all (5, A) G Zj for which there exists an elementary 

solution y of fl2]). This is a smooth irreducible algebraic curve in C 2 , [U|2]. In 
this paper we describe Zj (R) = Zj^flR 2 . We prove a result announced 
in [6]: 

Theorem 1. For n > 0, (R) consists of [n/2] + 1 disjoint analytic 

curves T nm , < m < [n/2] (analytic embeddings o/R to R 2 / 

For (6, A) G r njm , t/ie eigenfunction has n zeros, n — 2m of them real. 

If n is odd, then b — > +oo on both ends of each curve T n ^ m . If n is even, 
then the same holds form < n/2, but on the ends ofT nn /2 we have b — > ±oo. 

// (b, A) G r n m , (b, n) G r n>m+ i and b is sufficiently large, then n > A. 

This theorem establishes the main features of Z^f(R) which can be seen 
in the computer-generated figure in [3]. Similar results were proved in [5] for 
two other PT-symmetric eigenvalue problems. 

Our theorem parametrizes all polynomials P of degree 4 with the property 
that the differential equation y" + Py = has a solution with n zeros, n — 2m 
of them real [101 d IS] - 

Suppose that (b, A) G Zj ES (R). Then the corresponding eigenfunction y 
of ^2} can be always chosen real. Let y\ be a real solution of the differential 
equation in (J2]) normalized by y\{x) — > as x — > +oo, x G R. Then y x 
is linearly independent of y. Consider the meromorphic function / = yjy\. 
This function has no critical points in C, and the only singularities of f~ x are 
six logarithmic branch points. A meromorphic function in C with no critical 
points and whose inverse has finitely many logarithmic singularities is called 
a Nevanlinna function. All Nevanlinna functions / arise from differential 
equations y" + Py = 0, where P is a polynomial by the above construction: 
/ is a ratio of two linearly independent solutions of the differential equation. 

Consider the sectors 

Sj = {te ie : t > 0, \B- nj/3\ < tt/6}, j = 0, . . . , 5. 

The subscript j in Sj will be always understood as a residue modulo 6. 
Function / has asymptotic values oo, 0, c, 0, c, in the sectors S , . . . S 5 , where 
c G C. It is known that / must have at least 3 distinct asymptotic values [9], 
so c 7^ 0, oo. Function / is defined up to multiplication by a non-zero real 
number, so we can always assume that c = e llB , < < n, where the points 
and 7r can be identified. The asymptotic value c is called the Nevanlinna 



2 



parameter. There is a simple relation between c and the Stokes multipliers 

mi- 

The sectors Sj correspond to logarithmic singularities of the inverse func- 
tion Thus f^ 1 has 6 logarithmic singularities that lie over 4 points if 
c 7^ c, or over 3 points if c = c. 

The map (6, A) (->• /3 (mod7r), Zj ES (R) — > R is analytic and locally 
invertible [HI [2] , so /3 can serve as a local parameter on the real QES spec- 
tral locus. To obtain a global parametrization one needs suitable charts on 
Z® ES (R), where this map is injective. 

To recover /, one has to know the asymptotic value c and one more piece 
of information, a certain cell decomposition of the plane described below. 
Once / is known, b and A are found from the formula 

' = -^-^ + 2Jz-X). (3) 

Now we describe, following |4J, the cell decompositions needed to recover / 
from c. Suppose first that c ^ R. 




Fig. 1. Cell decompositions $ and $' of the sphere (solid lines). 



Consider the cell decomposition $ of the Riemann sphere C shown by solid 
lines in the left part of Fig. 1. It consists of one vertex at the point 2, 
three edges (loops j c , 7c and around non-zero asymptotic values) and 
four faces (cells of dimension 2). The faces are labeled by the asymptotic 
values 0, c, c, 00. Label is not shown in the picture. The face labeled is 
the unbounded region in the picture. (The point 1 in the figure is neither a 
label, nor a part of the cell decomposition. It will be needed, together with 
the dashed line L, for the limit at /3 = 0.) As 

/ : C \ f-\{0, 00, c, c}) -> C \ {0, 00, c, c} 
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is a covering map, the cell decomposition $ pulls back to a cell decomposition 
\I/ of the plane. 




Fig. 2. Two examples of the cell decomposition ^ of the plane (solid lines). 
Both eigenfunctions have two zeros, none of them real. 

Examples of \I/ are shown with solid lines in Figs. 2 and 4 (left). The faces 
of \I/ are labeled with the same labels as their images. Non-zero labels of 
bounded faces are omitted in the picture. The reader can restore them from 
the condition that labels around a vertex must be in the same cyclic order as 
in Fig. 1 (left, solid lines). The labeled cell decomposition \1/ defines / up to 
a pre-composition with an affine map of C. Two cell decompositions define 
the same / if they can be obtained from each other by a homeomorphism 
of the plane preserving orientation and the labels. Such cell decompositions 
are called equivalent. 

By replacing multiple edges of the 1-skeleton of \I> with single edges and re- 
moving the loops, we obtain a simpler cell decomposition T whose 1-skeleton 
is a tree, which we denote by the same letter T. The cell decomposition \I/ is 
uniquely recovered from its tree T embedded in the plane, [3]. The faces of T 
are asymptotic to the sectors Sj and the label of each face is the asymptotic 
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value in Sj. Two faces with a common edge cannot have the same label. 
The cell decomposition T is invariant under the reflection in the real axis, 
with simultaneous interchange of c and c. It is easy to classify all possible 
embedded planar trees T with labeled faces that satisfy these properties. 
They depend on two integer parameters k and I > 0. These trees form two 
families, X^i, k > 0, 1 > and X^i, k < 0, 1 > 0, shown in Fig. 3. Integers 
\k\ and I are the numbers of edges between ramification vertices, as shown 
in Fig. 3. 






Fig. 3. Trees X kjl . 

Cell decompositions in Fig. 2 (solid lines) correspond to the trees X 0j i and 
Xi i. Cell decomposition in the left part of Fig. 4 (solid lines) corresponds 
to the tree in the right part of Fig. 4. 

Parameters of the trees Xkj can be interpreted as follows: 

k~ := min{— k, 0} 

is the number of real zeros of /, and 21 is the number of non-real zeros. So 
the total number of zeros is n = 21 + k~. 

Functions / corresponding to the trees Xk t i, k > 0, have 21 zeros, none 
of them real. Zeros of the eigenfunction y coincide with those of /. 

For given n, the number of trees X^ with k < 0, 21 — k = n is 

(n + l)/2 when n is odd, and n/2 when n is even. (4) 



5 



Every tree Xk,i and every G (0, tt) defines a meromorphic function / satis- 
fying ([3]) with J = 21 + k~ + 1 and some (6, A) depending on (3, k and Z. This 
follows from a result of Nevanlinna [9], see also [I]. From this function /, 
the coordinates of a point (b, A) on the real QES spectral locus are recovered 
from the Schwarz equation (j3J). Thus we have a map F : (T, (3) i— >• (6, A) 
which we call the Nevanlinna map. This map is of highly transcendental 
nature: construction of / from T and (3 involves the uniformization theorem. 
We refer to [U [5J |9] for details. 

Each of the trees from our classification defines a chart of Z^ ES (IV). To 
obtain the global parametrization of Z® ES (R), we only have to find out how 
these charts are pasted together. We will see that the boundaries of our 
charts correspond to the values c = ±1. When c is real, we can use, instead 
of <£>, the cell decomposition $i of the sphere with two loops, and the 
loop L around c shown with the dashed line in the left part of Fig. 1. 

Proof of Theorem 1. We begin with the charts X^j, k < 0. We show that 
in these charts the limits as (3 — )■ 0, n do not belong to the spectral locus. 
This is proved by the arguments similar to those in [5J, Thm. 4.1]. 

Lemma 1. For k < and I > 0, the limit of the Nevanlinna map is 

lim F(X k>h f3) = oo0 

Proof. When — > 0, we have c — > 1, c — ^ 1. Suppose by contradic- 
tion that F(Xkj,(3) has a limit (bo, Ao). Then there is a limit function /o ; a 
solution of the Schwarz equation ([3]) with these parameters 60 and Ao- Mero- 
morphic function / has three asymptotic values, 0, 1, 00, and we are going to 
find the corresponding cell decomposition. Let $1 be the cell decomposition 
of the Riemann sphere with one vertex at the point 2 and two loops, 700 and 
L (see Fig. 1, left). Let ^ = 

It is easy to construct ^ from the original cell decomposition First, 
removing preimages of 7 C and 7c, we obtain the cell decomposition \lZoo, the 
preimage of the loop around 00 in $. It is shown with the bold solid lines in 
Fig. 4. 

Next, for each vertex v of ^ consider the path L v consisting of the edge 
of / _1 (7c) starting at v and ending at some vertex v', followed by the edge of 
/ _1 (7c) starting at v' and ending at some vertex v" . Then the edge of f^ 1 {L) 

^^Here 00 refers to a point added to the (6, A)-plane in the one-point compactification. 
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from v to v" is homotopic to L v in the complement of ty^. The new edges 
are shown with dashed lines in Fig. 4. The resulting cell decomposition is 
equivalent to 

Let V be the set of the vertices of ^ contained in the boundary of the 
sector S3. It is connected to the rest of the vertices of \& only at one vertex 
(v ' in Fig. 4, left) which is also at the boundary of both sectors S2 and S4. 
The dashed line replacing the edges of \P that connect v' to the two adjacent 
vertices of V (v and v" in Fig 4, left), goes from v to v" . All other dashed 
lines connect the vertices of V \ {v'} with the other vertices from the same 
set. Hence V \ {v'} is the set of vertices of a connected component of the 
1-skeleton of the cell decomposition This contradicts our assumption 
that = /o ($1), since the 1-skeleton of / _1 (<l ) i) must be connected. This 
contradiction proves the lemma. 




Remark. Consider all meromorphic functions with no critical points and 
at most 6 asymptotic values. These functions / are defined by their asymp- 
totic values and cell decompositions. Assume that one vertex t> of ^ is placed 
at z = and normalize so that f'{0) = 1. The class of normalized functions 
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obtained in this way is compact [12]. Let f v — > f be a converging sequence Jl 
The 1-skeletons of the corresponding cell decompositions ^(V) converge to 
the 1-skeleton of the cell decomposition ^(0) as embedded graphs with a 
marked vertex. If two asymptotic values collide in the limit, one has to use 
the procedure described in the proof of Lemma 1: replacing two loops by one 
loop. The limiting cell decomposition obtained in Lemma 3 suggests that 
the eigenvalue problem (T5]) tends to a harmonic oscillator when c — > 1, the 
fact we'll later prove by different arguments. 

Lemma 2. For k < and I > 0, 

lim F(X kh (3) = oo. 

Proof. When (3 — > n, we have c — > — 1, c — > — 1. Suppose by contra- 
diction that F(Xfc;,/3) has a limit (6 ,Ao). Then there is a limit function 
fo, a solution of the Schwarz equation Q with these parameters b and A . 
Meromorphic function fo has three asymptotic values, 0, —1, oo, and we are 
going to find the corresponding cell decomposition. 

To do this, it is convenient to choose another cell decomposition $' of 
the Riemann sphere, shown in the right part of Fig. 1 (solid lines). When 
c — > —1, $' collapses to where the two loops and 7^ are replaced with 
a single loop L' around —1 (dashed line in Fig. 1, right). 

We need the transition formula from \1/ = / _1 ($) to = /" 1 ($ / ). This 
formula is obtained by combining the two decompositions (see Fig. 5) and 
expressing the loops of $' in terms of the loops of $. 



2 Uniform convergence on compact subsets in the plane, with respect to the spherical 
metric in the target sphere. 
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Fig. 5. Two cell decompositions of Fig. 1 combined. 



The formulas, using notations in Fig. 5, are: 

7oo = a/3, i OQ = /3 a, y' c = Pjc/3' 1 , i s =a~ x 'y s a. (5) 

Here the product should be read left to right. Similar formulas were obtained 
in [5] in the proof of Theorem 4.1. Application of these transition formulas 
to the cell decomposition \1/ of type X_i 5 i is illustrated in Figs. 6,7. 




Fig. 6. Transition formulas (JSJ) applied to the cell decomposition in Fig. 4. 
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Fig. 7. Cell decomposition 



The same arguments as in Lemma 1 show that the 1-skeleton of the degen- 
eration of ty' as (5 — » 7r is not connected. This proves the lemma. 

Lemmas 1 and 2 show that for k < 0, the charts Xfc >m with 2m — k = n 
cover connected components of Z®+f(R), each parametrized by j3 £ (0, it). 
We call these components r n>m . These are simple disjoint analytically em- 
bedded curves in R 2 . 

When (3 — > 0, ir we must have 6 — >■ ±oo. We'll show below that b — > +oo 
on both ends of T n ^ m when > 0. 

When n is odd (that is J is even), these curves T n ^ m constitute the whole 
spectral locus Z^ H E 1 S (R). 

Now consider the part of the spectral locus covered by the charts Xk,i, k > 
0. This part is present only when n = 21 is even. 

Lemma 3. For k > and I > 0, we have 

lim F{X kh (5) = lim F(X k+hl , (5). (6) 

p— >7T ft— >0 

and 

KmF{X 0tl ,(3) =oo. (7) 
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Proof of Lemma 3. This is similar to the arguments in Lemmas 1 and 2. 
Computation is illustrated in Figs. 2, 8, 9 and 10. 




In the left part of Fig. 8 we use \I/ from Fig. 2, right. It corresponds to the tree 
Xi i in Fig. 8, right. In Fig. 9 the circles denote the vertices of (preimages 
of the vertex of $') and the dotted lines correspond to the preimages of j' c 
and 7^ determined from ([5]). The preimages of 7^ and 7^ coincide. They are 
shown with the bold solid line. Removing the preimages of 7 C and 7c (thin 
solid lines in Fig. 9) and the vertices of we obtain the cell decomposition 
shown in Fig. 10 (left) corresponding to the tree X 2 ,i (right). 
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Fig. 9. Passing from \l/ to 
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Thus, for k > 0, the cell decomposition of the plane obtained from 
Xkj in the limit /3 — > n as the preimage of is equivalent to the cell 
decomposition obtained from Xk+ij in the limit — > as the preimage 
of $i- Since = — $i, Nevanlinna theory implies that the corresponding 
functions / and /' satisfy /' = — /. Hence these two functions correspond to 
the same point of Zj (R). 

The proof of (|7j) is similar to that of Lemma 1. This completes the proof 
of the lemma. 

Now we continue the proof of Theorem 1. 

For even n = 21, charts Xk t i, k > parametrize segments of one curve in 
the real QES spectral locus, and we call this curve r„ jn / 2 . We parametrize the 
curve r„. n /2 by the real line, so that the number k decreases, thus the right 
end of r n n / 2 is parametrized by the chart X n / 2 . So when the parameter 
t G Ron r njn / 2 tends to +oo, the asymptotic value c = exp(z/3) tends to 
1. On the other hand, when t — > — oo on r n n / 2 the asymptotic value c does 



13 



not have a limit; it oscillates, passing each point of the unit circle infinitely 
many times. 

The curves r n>m are disjoint. Indeed, different cell decompositions give 
different functions /. This proves the first two statements of Theorem 1. 

Now we deal with the asymptotic behavior of our curves i\ m . We use the 
rescaling of fl2]) as in [B]. The QES spectral locus is defined by a polynomial 
equation Q n+1 (b,X) = which is of degree n + 1 in A. So on a ray b > 
bo there are n + 1 branches Xj(b). In [SJ Eq. (25)], we found that all Aj 
have asymptotics X(b) ~ b 2 + 0(y/b), b — > oo, and as b — > +oo, each QES 
eigenfunction yj tends to some eigenfunction Yg of the harmonic oscillator 

- Y" + 4z 2 Y = nY, Y(it) ->■ 0, t -> ±oo. (8) 

The eigenvalues of this harmonic oscillator are fie = 2(2£ + 1), £ = 0, 1, ... . 

Only one of the eigenfunctions yj can tend to a given Xe, and the corre- 
sponding eigenvalue satisfies 

Aj(6) = 6 2 + (/i £ -2J + o(l))A 6^+oo. 

It follows that all Aj are real. The graph of each Xj is a part of a curve T n m , 
and each T n m has only two ends. 

Now we consider the degeneration of the X 0j z chart with / > 0, the chart 
which parametrizes the right end of r n n / 2 , n = 21. On the left end of r n n / 2 , 
where t — > — oo in the parametrization described after Lemma 3, there are 
infinitely many points r n>n / 2 (£fc) which belong to the real QES locus, and 
where the asymptotic value c is real. It was proved in [6] that these are 
exactly those points where Zj ES (R) crosses the non-quasi-exactly solvable 
part of .Zj(R), and these points correspond to bk — > — oo. 

So only on one end of T n ^ n /2 (where t — >■ +oo) we can have b — > +oo. 
On the other hand, each r n>m , m < n/2 contains at most two graphs of Aj. 
According to (j4j), the total number of these graphs Xj is n+ 1, and the total 
number of curves T n ^ m is (n + l)/2 when n is odd, and n/2 + 1 when n is 
even. It follows that, when n is odd, each T n ^ m contains two graphs of Aj. 
When n is even, each r nm except one contains two graphs of Aj, while the 
exceptional component r n n / 2 contains one graph of Aj. 

Thus b — > +oo as c — > ±1 in the X/^-charts with k < 0, which proves the 
third statement of Theorem 1. To prove the last statement, we study zeros 
of the eigenfunctions as b — > +oo. 
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The eigenfunction Y e of flH]) corresponding to the eigenvalue /i£ has exactly 
I zeros on iR and no other zeros in C. One of these zeros is real iff I is odd. 

The trees corresponding to Yi are constructed similarly to those corre- 
sponding to y, using the two loop cell decomposition of the sphere, consisting 
of 7oo and the dashed loop in Fig. 1, left. 




Fig. 11. The tree corresponding to Y 3 . 

For general results on convergence of Nevanlinna functions like our / we refer 
to [LJ. 

When b — > +00, each QES eigenvalue X(b) must tend to some /i^, and 
the corresponding QES eigenfunction tends to Yg. Suppose that A (6) G T n ^ m 
with m < n/2. Then the tree corresponding to X(b) is X/,^ m , k < 0, 2m — 
k = n. From the arguments in the proofs of Lemmas 1 and 2 (see Figs. 4, 
7), degeneration of the cell decomposition \1/ corresponding to such a tree 
has a connected component with 2m bounded faces when /3 — ?• and with 
2m + 1 bounded faces when /3 — > tt. This implies that the corresponding 
eigenfunction can only converge to Y 2m as (3 — > and to Y 2m+ i as — > tt. 

If n is odd, these curves constitute the whole QES locus. If n is even, 
there is one more branch A (6) of the QES locus for large positive b, the right 
end of corresponding to the tree X 0in /2. From the proof of Lemma 3 

(see Fig. 2, left) degeneration of the cell decomposition \l/ corresponding to 
such a tree has a connected component with n bounded faces when /3 — > 0. 
This implies that the corresponding eigenfunction can only converge to Y n . 
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So the ordering of the ends of the curves r„ m corresponds to the natu- 
ral ordering of the first n + 1 eigenvalues of the harmonic oscillator. This 
completes the proof. 
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